Proceedings of the 363. WE-Hcraeus Seminar on: "Neutron Stars and Pulsars" (Posters and contributed talks) 
Physikzentrum Bad Honnef, Germany, May. 14 - 19, 2006, eds. W. Becker, H.H. Huang, MPE Report 291, pp. 36- 39 



Gamma-ray Pulsar Simulations for the Gamma-ray Large 
Area Space Telescope (GLAST) 



M. Razzano 1 2 



o 
o 

o 

Q 

(N 



> 

m 

(N 
O 

a: 

6 

CO 

C3 



1 Istituto Nazionale di Fisica Nucleare, Pisa sect., Largo B.Pontecorvo 3, 56100 Pisa, Italy, emailmassimiliano. razzano@pi.infn.it 

2 on behalf of the GLAST LAT Collaboration 



Abstract. We present here the current status of pulsar 
simulations for the Large Area Telescope, the main instru- 
ment on the GLAST mission. We present PulsarSpectrum, 
a pulsar simulator that can reproduce with high detail 
gamma-ray emission from pulsars. PulsarSpectrum takes 
into account advanced timing effects, e.g. period changes 
with time, barycentering effects and glitches. Other ancil- 
lary tools have been built to provide the simulator with a 
realistic population of pulsars and their ephemerides. All 
these tools are currently used in the GLAST collabora- 
tion for testing the LAT Science Analysis Tools and for 
studying LAT capabilities for pulsar science. They have 
been used also for generating a simulated pulsar popu- 
lation for the Data Challenge 2 (DC2), one of the most 
important milestones in the development of the GLAST 
software. During the DC2, scientists analyzed a set of 55 
days of simulated data in order to validate LAT Monte- 
Carlo, study instrument response functions, exercise anal- 
ysis tools and study LAT capabilities. This contribution 
also contains results of some analysis performed during 
DC2 on EGRET pulsars. 



1. Introduction 

Pulsars arc among the most exciting gamma-ray sources in 
the Universe and can serve as unique sites for the study of 
emission processes in extreme physical environments. The 
Gamma ray Large Area Space Telescope (GLAST) will in- 
crease dramatically our knowledge of gamma ray pulsars 
physics. In particular the Large Area Telescope (LAT), 
the main GLAST instrument, will provide more detailed 
observations of the known gamma-ray pulsars and poten- 
tially will discover many new pulsars that emit gamma 
rays. To better understand the capabilities of GLAST for 
pulsar science we developed PulsarSpectrum, a program 
that simulates gamma ray emission from pulsars and a set 
of ancillary tools that can create simulated pulsar data to 
be used with PulsarSpectrum. This simulator can be eas- 
ily interfaced with the MonteCarlo software that simulates 
the response of the LAT. 

Today seven high-confidence and three low-confidence 




Fig. 1. The GLAST Large Area Telescope. An incident 
gamma ray converts into a e~e + pair 



gamma-ray pulsars are known and this numbers will in- 
crease with the observations of GLAST, scheduled for 
launch in the fall of 2007. 

The GLAST Large Area Telescope (LAT) (Fig. Q}, is 
a pair conversion telescope based on the most advanced 
high energy detectors. It consists in a precision silicon 
strip tracker, an hodoscopic calorimeter and a segmented 
anticoincidence shield for particle background rejection. 
The LAT high sensitivity (2xl0~ 9 ph/cm 2 /s) and effec- 
tive area (>8000 cm 2 ) will permit the discovery of a lot of 
new pulsars: the estimates range between tens to hundreds 
depending upon the theoretical scenario adopted. More- 
over the low dead time of the detector (20 fj,s) will allow 
the detailed reconstruction of pulsar lightcurves. One of 
the most exciting possibilities of the LAT will be the cov- 
erage of the energy window from 30 GeV up to 300 GeV, 
a still unexplored range. At these energies the theoreti- 
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Table 1. The PulsarSpectrum parameters. 



ran 



Type of Parameter 



Parameter 



Observational" : 
Model-dependent : 



RA, Dec, flux, 
t ,f(t ),f'(t ),f"(t ), 
Model, E 

Random seed, 
5 free parameters, 



to indicates the epoch of the ephemerides 



cal models make different previsions on the high energy 
spectral cutoff and the spectral coverage of LAT will be of 
primary importance for constraining and discriminating 
among the models. In order to study the LAT response to 
specific gamma ray sources, various simulation packages 
have been developed. 

2. The PulsarSpectrum simulator 

2.1. General overview 

The basic idea of PulsarSpectrum is to construct a bidi- 
mensional histogram representing the differential flux vs. 
energy and pulsar phase. This histogram contains all the 
basic informations about lightcurve and spectrum. How 
it is built depends upon the model we want to use: for 
example a phenomcnological model, based only on obser- 
vations, or a physical one, based on a specific theoretical 
model. Currently two models have been included. The first 
model implemented was phenomcnological, since it is more 
flexible and completely independent from the theoretical 
emission scenario. A second model recently implemented 
allow the user to simulate pulsars with an arbitrary pho- 
ton distribution in phase and in energy. This model can 
in fact use a 2D histogram that can be provided by other 
external programs. 

The input parameters of the simulator can be divided in 
two categories and are listed in the Table [TJ 

— Observational parameters, which characterize the gen- 
eral parameters of the pulsar; 

— Model- dependent parameters, that define which model 
will be used for simulation and the set of parameters 
used by this model. There are 5 free parameters that 
leave to the user the freedom to adjust the model. 

All parameters are placed in two specific data files used 
by both PulsarSpectrum and the LAT simulation tools. 
PulsarSpectrum creates the lightcurve and the spectrum 
from these parameters and combines them to obtain 
a two-dimensional matrix that represents the flux in 
ph/m 2 /s/keV. An example of such an histogram for a sim- 
ulated pulsar similar to Vela is in Fig. [2J The photons are 
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Fig. 2. An example of a 2D histogram generated by Pul- 
sarSpectrum for simulating a pulsar with a spectrum sim- 
ilar to Vela pulsar 



then extracted such that the interval between two subse- 
quent photons is determined by the flux integrated over 
the energy range of interest. 

2.2. The phenomenological model 

The phenomenological model allows the user to generate 
pulsar lightcurves in a general way using a single or dou- 
ble Lorenzian peak profile whose shape is determined from 
random generated numbers. The lightcurve can be gen- 
erated alternatively from a user-provided profile. This is 
useful for simulating the already observed gamma ray pul- 
sars. The spectral shape is assumed to be a power law with 
exponential cutoff (according to [TJ), as in the observed 
gamma-ray pulsars and can then be modeled as: 



dN r-, E xa 



(1) 



The normalisation constant K is determined by the pho- 
ton flux in the range 100 MeV- 30 GcV, in order to have 
flux compatible with the fluxes in the 3 rd EGRET Cata- 
log. The other spectral parameters can be varied; the val- 
ues for the EGRET pulsars are obtained from real data 
by fitting procedures (See e.g. [I]). As an example, in 
Fig. [2] is presented an histogram for a Vela-like pulsars, 
i.e. a simulated pulsar with characteristics similar to the 
Vela pulsar. The lightcurve has been derived from EGRET 
observations while the spectral parameters from a multi- 
wavelength fit as reported in [TJ are displayed in Fig. [3] 



2.3. Timing issues 

Once the differential flux histogram is created the time 
interval between two subsequent photons is computed ac- 
cording to the flux. The strategy adopted is to compute 
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Fig. 3. The input spectrum for a Vela-like pulsar using 
phenomcnological model. The spectral parameters have 
been obtained from pQ. The crosses indicate the distribu- 
tion of the extracted photons. 




Fig. 4. The count map of the simulatd data of the LAT 
DC2. The green circles indicate the EGRET pulsars, part 
of the whole pulsar population simulated using Pulsar- 
Spectrum 

the position of the Sun. The accuracy for the computation 
of these corrections is hard-coded in the program. 



the mean photon rate and then to calculate the inter- 
val to the next photon according to Poisson statistics. 
Then a secondary correction At is summed to this in- 
terval in order that the photon distribution follow the in- 
put lightcurve. The interval between two photons is com- 
puted assuming that the pulsar period does not change 
with time and the photons arrival times are computed 
into a reference system fixed relative to stars. This is not 
the " real world" . Pulsar timing is affected by more com- 
plicate effects, as (1)- The motion of GLAST through the 
Solar System and the relativistic effects due to gravita- 
tional well of the Sun (see l2.3.Ijk (2)- Period changes with 
time (see I2.3.2[) . For pulsars in binary systems an addi- 
tional modulation to the orbital period should be taken 
into account. For a precise pulsar simulator intent to pro- 
duce a realistic list of photon arrival times we need to 
include all these effects (to transform to the observational 
frame). All these procedures arc now implemented in the 
code and only the binary demodulation is not yet imple- 
mented. The real arrival time of a photon from a pulsar 
must be first barycentered and then phase assigned. 

2.3.1. Barycentric effects 

The first step to analyze pulsar data is the conversion 
from the arrival times at the spacecraft, usually expressed 
in Terrestrial Time TT or TAI, to the arrival times at 
the Solar System barycenter, expressed in Barycentric Dy- 
namical Time TDB. Taking into account both the mo- 
tion of spacecraft through space and the general relativis- 
tic effects due to the gravitational field of the Sun (i.e. 
Shapiro delay), the simulator computes the opposite of 
the barycentric correction by considering the position of 
the Earth and of the spacecraft in the Solar System, and 



2.3.2. Period change and ephemerides 

The rotational energy of a radio pulsar decreases with time 
and hence the period increases with time. For gamma-ray 
pulsar science the radio ephemerides arc fundamental for 
assigning the correct phase to each photon. If we know 
the frequency f(tg) and its derivatives /(to) and /(to) at 
a certain time to, known as epoch, the phase is then: 



<f>(t)=int[f{t ){t-to)- 



|/(*o)(*- 



1 



toy+-f(t ){t-t y 

o 



(2) 

The interval between two photons must be also corrected 
for this effect. In the parameters file the user can specify 
a set of ephemerides with the relative epoch of validity 
expressed in Modified Julian Date. The simulator then 
computes the opportune arrival time such that, after ap- 
plying the barycentric corrections and then the Eq. [2j the 
correct phase is obtained. 

3. Gamma-ray pulsars in the LAT Data 
Challenge 2 

Pulsar Spectrum have been used by the GLAST LAT Col- 
laboration for simulating pulsars in the LAT Data Chal- 
lenge 2. The Data Challenge 2 (DC2) was an important 
milestone for the preparation of the analysis and simu- 
lation software for the GLAST mission. For the DC2 a 
simulated 55-days long LAT observation of the whole sky 
have been generated, and the simulated data have been 
provided to the scientists of the Collaboration. The main 
goals of the DC2 were to study and validate the LAT Mon- 
teCarlo, to exercize the Analysis Tools under development 
for GLAST and to better study the LAT Instrument Re- 
sponse Functions. Using DC2 simulated data, a group of 
LAT scientists have worked on these from the begin of 
March to the end of May (date of official DC2 Closeout 
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Fig. 5. The phase distribution of the photons above 100 
McV from the Vela pulsar using DC2 simulated data. 
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Fig. 6. The phase distribution of the photons above 100 
MeV from the Geminga pulsar using DC2 simulated data. 



meeting) and beyond. 

One of the most important efforts in the DC2 was the com- 
position of an high-detailed simulation of the gamma-ray 
sky (see Fig. HJ) - This similated sky was built by using a set 
of programs capable to simulate the different gamma-ray 
sources. Not only the current known gamma-ray sources 
discovered by EGRET have been included, but also a lot of 
new sources, in order to push down the sensitivity limit of 
the detected sources. Also possible new classes of gamma- 
ray sources have been included, such X-ray binaries and 
microquasars. 

The pulsar component was entirely simulated by Pulsar- 
Spectrum and was composed by some subclasses of pul- 
sars. There are the six gamma-ray pulsars detected at 
EGRET energies, displayed in Fig. HJ and some other 
real radio-pulsars that are within error box of some 3EG 
sources. The remaining pulsars were not actually real pul- 
sars, but simulated pulsars that have characteristics sim- 
ilar to the current known pulsars. These were divided in 
two subgroups, a set of isolated pulsars and a set of mil- 
lisecond pulsars. A database containing the timing solu- 
tion for some pulsars has been provided to DC2 users. 
However for some gamma-ray pulsars the timing solution 
have not been provided, in order to leave to the DC2 peo- 
ple the possibility to study tccniqucs for finding Geminga- 
like pulsars with no radio counterparts. In Fig. [5] and [5] 
are displayed the recontructed phase curves of the simu- 
lated Vela and Geminga pulsar using DC2 simulated data. 

These curves showed a profile very similar to the real 
EGRET profile, but this is not a suprise since the EGRET 
data have been used for shaping the input lightcurve for 
EGRET pulsars present in the DC2. With comparison 
of EGRET it is remarkable the better statistics achieved 
with the LAT in 55 days in comparison to the EGRET 
curves, as reported e.g. in [2] 



4. Conclusions 

Pulsar simulations are very important to study the 
GLAST LAT capabilities to probe these gamma-ray 
sources. We present here the latest results in pulsar simu- 
lations using the PulsarSpectrum simulator, able to repro- 
duce timing and spectral features of gamma-ray pulsars 
with high detail. PulsarSpectrum has been successfully 
used by the GLAST collaboration for testing the function- 
ality of the LAT Science Analysis Tools, a set of analysis 
programs specifically designed to analyse the LAT data af- 
ter launch. An important event in the preparation of the 
GLAST mission was the Data Challenge 2, during what 
a 55-days LAT observation of the full sky have been gen- 
erated. Pulsars included in DC2 sky have been simulated 
with this simulator, providing LAT scientists simulated 
data from a realistic gamma-ray pulsar population. The 
DC2 has been also a good chance to exercize the pulsar 
simulator and to show that PulsarSpectrum can be used 
by the GLAST Collaboration for better exploring the LAT 
detecting capability on pulsars. 
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